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Abstract 



i-C | We consider pair production and decay of fundamental unstable particles in 

the framework of a modified perturbation theory (MPT), which treats resonant 
^■i contributions of unstable particles in the sense of distributions. The cross- 

r- 1 [ sections for the top-quark pair production and for the VF-boson pair production 

in e + e~ annihilation are calculated within the NNLO in models that admit 
exact solutions. In both cases an excellent convergence of the MPT is detected 
at the energies close to and above the maximum of the cross section. In the case 
" of W^-boson pair production a precision of the description at the ILC energies 

■ is ensured at the level of one per-mille or higher. 

The processes of pair production and decays of unstable fundamental particles, 
such as top quarks and W bosons, play an important role for testing the Standard 
Model and for searching physics beyond. In the case of colliders subsequent to LHC a 
description of such processes must be made generally with the NNLO accuracy. This 
implies that not only the gauge cancellations and unitarity should be maintained, 
but also suitably high accuracy of computation of resonant contributions must be 
provided. Unfortunately, the existing methods can provide only the NLO precision 
of the description of the cross- sect ion. This is the case with the double pole approx- 
imation (DPA) successfully applied at LEP2 [lj or with the complex-mass scheme 
(CMS) [2] intended mainly to ILC |3J. The pinch-technique method, another for a 
long time developed approach, in principle can provide the NNLO precision, but to 
maintain the gauge cancellations it requires a huge volume of calculations of extra 
contributions that pertain formally to the next level of the precision, which is im- 
practical [5]. So alternative approaches are required for systematic calculations at 
the NNLO. A modified perturbation theory (MPT) [6], [7J E] is a suitable approach 
for solving this problem. Its main feature is the direct expansion of the probability 
instead of amplitude in powers of the coupling constant with the aid of distribution- 
theory methods. As the expansion is made in powers of the coupling constant and the 
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object to be expanded is gauge invariant, the gauge cancellations in the MPT must 
be automatically maintained. Nevertheless, the accuracy of the description in MPT 
requires an examination. In order to do that, numerical simulations are necessary in 
the framework of MPT. 

Since in the case of pair production of unstable particles the most crucial are the 
double-resonant contributions, we consider initially only these contributions in the 
framework of the model simulations. (Generally, the single-resonant contributions 
may be considered, as well, this is not a problem in the MPT approach |8J.) In the 
case of e + e~ annihilation, the corresponding total cross-section has the form of a con- 
volution of hard-scattering cross-section with the flux function. The hard-scattering 
cross-section is an integral over the virtualities of unstable particles of exclusive cross- 
section multiplied by factor standing for soft massless-particles contributions, 

a(s) = [ — 4>(s'/s;s) a(s') , a(s') = ff dsi ds 2 a exd (s', s x , s 2 ) (l+5 so/t ) . (1) 

J s min s J J 

The exclusive cross-section a exd is written as a product of Breit-Wigner factors p(si), 
some kinematic factors, and a function $, which is the rest of the amplitude squared, 



&excl(s,S 1 ,S 2 ) = 0^/s-y/sl - y/s^)\ / \{s,S 1 ,S 2 ) $(s; Si, S 2 ) pOl) p(s 2 ) • (2) 

Generally $ corresponds to one-particle irreducible contributions, and so it does not 
have singularities on the mass-shell of unstable particles. On the contrary, the kine- 
matic factors, which include the theta-function and the square root of the kinematic 
function A, have singularities. The BW factors if to naively expand them in powers 
of the coupling constant a generate non-integrable singularities, and this makes up a 
great problem because integrals in ([I]) become senseless. 

However, the singularities become integrable if to expand the BW factors in the 
sense of distributions. In this case the expansion of a separately taken BW factor is 
beginning with the 5-function which corresponds to the narrow-width approximation. 
The contributions of the naive Taylor expansion are supplied with the principal-value 
prescription for the poles. The nontrivial contributions are the delta-function and its 
derivatives with coefficients c n , which are polynomials in a that are determined by 
the self-energy of the unstable particle [6]. Within the NNLO, the expansion is as 
follows: 

p( 8 ) = MH\ s -M 2 + i:(s)\- 2 = 5(s-M 2 ) (3) 
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Here M is the renormalized mass, T is the Born width, S(s) is the self-energy of the 
unstable particle. Coefficients c n within the NNLO include 3-loop self-energy contri- 
butions and their derivatives determined on-shell. The structure of the contributions 
is such that in the OMS-type schemes of the UV renormalization the real self-energy 
contributions enter into the coefficients either without the derivatives or with the 
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first derivative only. This means that the relevant real self-energy contributions are 
determined by the renormalization conditions. In the case of unstable particles it is 
reasonable to use the OMS or pole scheme [HI EE] of the UV renormalization, whose 
inherent property is that the renormalized mass of unstable particle by definition co- 
incides with the real part of the pole of unstable-particle propagator, which is gauge 
invariant and scheme-independent. The coefficients c n in this scheme are determined 
as follows El: 



^2 . 2 

Co = — a — + a 
h 



ci = 0, c 2 = - a 2 If . (4) 



Here I k = ImS fc (M 2 ), E = a Ei+a 2 E 2 +ct 3 E 3 , and I[ = ImE^M 2 ). Simultaneously 
in the OMS scheme the ImE(M 2 ) coincides with the imaginary part of the pole of 
the propagator. This allows one to determine Ik order-by-order over the width of 
the unstable particle via the unitarity relations <xl\ = MTq, a 2 I 2 = MaT±, and 
a 3 I 3 = Ma 2 T 2 + T 3 /(8M) 

Unfortunately, expansion ([3]) has sense only if the weight in the integral is a regular 
enough function. In our case, however, the kinematic factors are not regular, which 
leads to a divergence in integrals (CQ) after the substitution of the expansions. At 
first glance, this brings up a question about the applicability of expansion (jHJ). Nev- 
ertheless, the kinematic factors may be analytically regularized via the substitution 
[A(s, S\, s 2 )] 1 ^ 2 — >■ [A(s, Si, s 2 )] u . With large enough v this imparts enough smooth- 
ness to the weight, and the singular integrals become integrable. Fortunately, after 
the analytic calculation of integrals the regularization may be removed without the 
loss of finiteness of outcomes. Moreover, the expansion remains asymptotic |8J. This 
completely salvages the applicability of the approach. 

The scheme of the analytical calculations is as follows. At first we proceed to 
dimensionless energy variables x, Xj (i = 1,2) counted off from thresholds, \fs = 
2M(1 + x/4), y/sl = M(l + X{/2). The hard-scattering cross-section then takes the 
form 

a(x) = // dxi dx 2 (x — Xi — x 2 )+ ; Xi, x 2 ) p(xi)p(x 2 ) • (5) 



Here (x — x\— x 2 )^_ = 6(x — x± — x 2 )(x — x± — x 2 ) u and tilde marks the dimensionless 
functions (factor l+5 so f t is included into the definition of $). Further, we substitute 
asymptotic expansions for p(xj), and consider at each (i — 1, 2) the contributions of 
PVx~ Hi and 5 <yril ~ l \x i ). Simultaneously, in each case, we represent the test function 
in the form of a double Taylor expansion over Xj truncated at the contributions of 
x - ni 1 , with a remainder, 



711-1712-1 kl k 2 

$(z;x ls x 2 ) = EEAn $ (fcl ' fe2) (x;0,0) + A$(x;x 1? x 2 ). (6) 

The higher powers of x, in the Taylor expansion will zero the 5^ ni ~ l \xi) and cancel 
the PVx~ ni . The remainder A$ is determined as the difference between $ and the 
Taylor expansion. In fact A$ is to be further expanded with respect to separately 
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x\ and X2, but for brevity we do not consider here this procedure explicitly (see 
details in [8]). Let us mention only that the final remainder produces a regular 
contribution to the integrand in formula (jSJ) , and the integrals of it can be numerically 
calculated at v = 1/2. At the same time, the contributions of Taylor are singular. 
However, the integral (jSJ) of them may be analytically calculated owing to simple 
(power) dependence on X{. After making the calculation and after putting v — > 1/2, 
the result appears in the form of a sum of regular and singular contributions with 
singular contributions being products of regular factors and power distributions of the 
type x 5 _/ 2 ~ n with integer n. It should be emphasized that at this stage of calculations 
the test function $ is determined by means of conventional perturbation theory, but 
the analytic calculations are made independently of the particular form of $. The 
convolution integral of the result can be numerically calculated. In particular, the 
integral of singular power distributions can be calculated by means of the formula 



where is a weight and N is a positive integer such that — N— 1 < Rez/ < — N. 

For carrying out further numerical calculations a double-precision FORTRAN 
code is written. The computation of regular integrals in this code is realized by 
Simpson method. Numerous indeterminate forms of the type 0/0 that emerge in the 
integrand due to the difference structures are resolved through the introduction of 
linear patches. The patches diminish the errors (numerical instabilities) that arise 
because of the loss of decimals near the indeterminacy points. At the same time, 
the errors generated by patches may be numerically estimated. Ultimately, the total 
errors caused by indeterminacies may be estimated, as well. The crucial point is 
that the errors because of the patches are increasing with increasing the sizes of the 
patches, while the errors because of the loss of decimals are decreasing. So there 
should be an optimum size of the patches when the sum of the errors is minimized. 
The minimization point must possess extremum properties, so that the result of the 
computation at this point must be stable with respect to varying the sizes of the 
patches. Furthermore, at the extremum point the sums of the errors of different 
kinds must be approximately equal each other (up to a coefficient of order one). So 
the order of the total error may be estimated by the order of the sum of the errors 
because of the patches. More details of how to do estimation of the errors, is found in 
[TT] . Eventually the adjusted estimate of the relative error of the computation of the 
NNLO approximation turns out to be less than 10~ 3 or 10~ 4 in two cases considered 
below. 

The physical models underlying the calculations are related to processes e + e~ — > 
7, Z ->- tt -> W + bW~b and e+e" ->• 7, Z ->- W + W~ -> 4/. Recall that our 
aim is to verify whether the MPT calculations are realizable and then to test the 
convergence properties of the MPT expansion. Having that in mind, we consider 
the test function <3> in both cases in the Born approximation. However, the self- 
energies in the denominators of propagators of unstable particles, we consider in 3-loop 
approximation (see details in [El [13]). This means that we can immediately check by 
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Figure 1: The total cross-section in the model for tt production and decay. 
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Figure 2: The total cross-section in the model for W + W production and decay. 



our calculations the convergence properties of the MPT expansion of the products of 
BW factors. Actually, this is sufficient for our purposes because the insertion of the 
loop corrections to the test function may be considered as the replacement of the test 
function by another one with additional factors a, a 2 , etc. Meanwhile, as the existing 
experience shows, the convergence properties of the MPT expansion are very weakly 
depend on the test function, but depend on the values of the corrections to the widths 
of unstable particles [T2J [T3] . As concerns the soft massless-particles contributions, 
we consider among them only universal ones. They are collected in the flux function 
and in the Coulomb factor. The flux function, we take into consideration in the 
leading- log approximation. The Coulomb factor, we consider in the one-gluon/photon 
approximation with specific resummation [2] that does not affect the BW factors. 
Notice that although the multi-gluon contributions generally are important in the 
case of the top quarks, we believe that at distance from the threshold a qualitative 
picture may be simulated in the one-gluon approximation. 

The outcomes of computations are presented in the Figures and in the Tables. In 
Fig. [TJ in the panels (a) the thick curve shows the behavior of the total cross-section 
in the model in the case of the top-quark pair production. The dotted, dashed, and 
continuous thin curves show the results of the MPT computations in the LO, NLO, 
and NNLO approximations, respectively. It is worth noting that the NNLO result 
almost coincides with the exact result in the region near and above the maximum of 
the cross-section. The distinction is visible in the panel (b) where the percentages 
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Table 1: tt production: the total cross-section in pb and 
in the model. 


in % with respect to exact result 


yfi (GeV) 
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O'NLO 


&NNLO 
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0.6724 


0.5687 


0.6344 


0.6698 




100% 


84.6% 


94.3% 


99.6% 
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0.2255 


0.1821 


0.2124 


0.2240 
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80.8% 


94.2% 


99.3% 


1500 
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Table 2: W + W production: the total cross-section in pb and in % with respect to exact 
result in the model. 


yfi (GeV) 


0" EXACT 
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100% 
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2.8286 
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2.8285 




100% 
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100.00% 


3000 


0.61023 


0.55733 


0.60625 


0.61026 




100% 


91.33% 


99.35% 


100.00% 
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with respect to the exact result are presented. In Fig. [2] the appropriate results are 
presented in the case of W pair production. The main difference at Fig. [2] is that 
already the NLO curve almost coincides with the exact result with the distinction 
visible in the panel (b) only. In Fig. [2] the appropriate results are presented in the 
case of W pair production. The main difference at Fig. [2] is that already the NLO 
curve almost coincides with the exact result with the distinction visible in the panel 
(b) only. In Fig. [3] the results for NLO and NNLO of Fig. [2jb are repeated with 
greater scale on vertical axis. In Tables [T] and [2] the outcomes are represented in the 
numerical form at the characteristic energies at ILC [3J. Note that the relative errors 
of the computation are less than 10~ 3 and 10 -4 in the cases of the top quarks and W 
bosons, respectively. Therefore the errors are omitted in the data of Tables. 

In conclusion, first of all we emphasize that the above results show in practice 
the existence of the MPT expansion in the case of pair production and decay of 
fundamental unstable particles. Secondly, the NLO and NNLO approximations in 
the MPT have very stable behavior at the energies near and above the maximum of 
the cross-section. The latter result to a large extent is model-independent since it 
is established in different models. (The latter point is discussed in more details in 
[T2j IT5] . Notice also that at lower energies, in particular near threshold, the mode of 
MPT that has been considered here becomes inapplicable, but there is another mode 
for MPT |8J.) Thirdly, at the ILC energies the NNLO approximation in the MPT give 
highly satisfactory results in numerical sense. Namely, in the case of the top-quark 
pair production it gives approximately a half-percent accuracy of the description of 
the cross-section, and in the case of W bosons does a per-mille accuracy. In fact, 
this is what is needed at the ILC. So, we conclude that MPT is a good candidate for 
support at the ILC the pair production and decay of fundamental unstable particles. 
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